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Abstract 

Despite the widespread presence of magnetic fields, their origin, evolution and role are 
still not well understood. Primordial magnetism sounds appealing but is not problem free. 
The magnetic implications for the large-scale structure of the universe still remain an open 
issue. This paper outlines the advantages and shortcomings of early-time magnetogenesis 
and the typical role of B-fields in linear structure-formation scenarios. 

1 Introduction 

Magnetic fields are everywhere in the universe pp. The Milky Way and many other galaxies 
possess fields of /uGauss strength, while ordered .B-fieids in the intracluster space suggest that 
structure formation could have been influenced by magnetic forces. In addition, Faraday rotation 
measurements at high redshifts indicate dynamically significant £>-fields in protogalactic clouds. 
Despite its widespread presence, however, the origin of cosmic magnetism remains a mystery 
and a subject of debate. The alignment of the galactic fields supports the dynamo-amplification 
idea, but dynamos require an initial seed field to operate [2]. These seeds could be the result 
recent (post-recombination) physics, or have a primordial (pre-recombination) originQ Deciding 
that is difficult, since the galactic fields have lost memory of their earlier history. In contrast, 
possible 5-fields in the intercluster space, or magnetic imprints in the CMB should provide 
much better insight. 

Primordial magnetism is attractive because it could explain all the magnetic fields seen in 
the universe, especially those in high redshift systems. Early magnetogenesis is not problem- 
free however pQ. Magnetic fields generated in the radiation era, namely between inflation and 
recombination, have too small coherence lengths and will distabilise the dynamo. Inflation 
can provide large-scale correlations, but i?-fields that survived a epoch of de Sitter expansion 
are typically too weak to sustain the dynamo. Primordial magnetic fields must also comply 
with certain constraints. These come from nucleosynthesis and, mainly, the CMB, which limits 
the current magnetic strength between 10 -6 Gauss (for random fields) and 10~ 9 Gauss (for 
homogeneous ones). Cosmological S-fields of /iGauus strength can affect structure formation, 
through their Lorentz force and the anisotropy that they introduce [5]. The Lorentz force, in 
particular, generates and affects all types of density inhomogeneities. The anisotropic nature of 
the field, on the other hand, makes it a source of shear and gravitational waves. 



1 At recombination, roughly 10 s years after the Big Bang, matter ceases being ionised and the universe becomes 
'transparent' to photons. Relic of that epoch is the Cosmic Microwave Background (CMB) radiation [3]. 
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2 Cosmological electromagnetic fields 



Consider a general spacetime with metric g a b, of signature (—,+,+,+), allow for a family of 
fundamental observers with 4- velocity u a (such that u a u a = — 1). The tensor h a b = g a b + u a u b 
projects orthogonal to u a into the observers' 3-D rest space. Together, h a b and u a introduce 
an 1+3 'threading' of the spacetime into time and space, while they decompose all physical 
and geometrical quantities into their timelike and spacelike parts. For instance, splitting the 
gradient of the 4- velocity, 

V&u a = - Qh ab + (J ab + iO a b ~ A a u a , (1) 

gives the kinematic variables [6j. These are the average volume expansion/contraction O = 
V a u a , the shear a ab = D( b u a }, the vorticity u> a b = D[ b u a ] and the 4-acceleration A a = ii a = 
u b VbU a (V a is the covariant derivative operator and D a = h a b Vb is its 3-D counterpart) H 
Similarly, the energy-momentum tensor of a general imperfect fluid, decomposes as [U] 

T ir } = P U a u b + PKb + 2<7( a U fe) + Tr ab , (2) 

where p, p, q a and Tt a b are respectively the density, the isotropic pressure, the energy flux and the 
anisotropic pressure of the matter. Also, q a u a = = 7r a bU b , with 7r ab = 7r ba and ir a a = 0. For a 
perfect fluid q a = = ir a b and for a barotropic medium p = p{p). Note that we use geometrised 
units, with c = 1 = 8ttG, throughout this paper. 

Relative to the fundamental observers, the electromagnetic field splits into an electric and 
a magnetic component, represented by the vectors E a and B a respectively (with E a u a = = 
B a u a ). Then, the stress-energy tensor of the Maxwell field reads as [5] 

<em) _I (t& , R 2A ,„ , I (-pi , R 2 



T ab ' = g \ E + B ) U * Ub + q{ E+ B J Kb + 2Q ^ Ub ) + Vab ' ^) 

with E 2 = E a E a , B 2 = B a B a , Q a = e abc E b B c being the Poynting vector and V ab = -E {a E b) - 
B( a B b ^ the anisotropic electromagnetic pressure. By construction Q a u a = = V a bU a , with 
V ab = V ba and V a a = 00 In other words, the Maxwell field corresponds to an imperfect fluid 
with p^) = (E 2 + B 2 )/2, p(^ = (E 2 + B 2 )/6, qa em) = Q a and ^ = 

The evolution of the electromagnetic field is monitored by Maxwell's equations. In the 
u a -frame, these decompose into a set of two propagation equations [5 J , 



E {a ) = 



<y a b + w«t - \ Ohab^j E b + e a b c A h B c + curLB a - J a , (4) 



2 - . ^b 



B {a ) = ( v a b + Uab ~ 3 <S>h ab ) B b - E a b c A b E c - curl^ a , (5) 



2 Round brackets denote symmetrisation and square ones antisymmetrisation. Angled brackets indicate the 
symmetric and trace-free part of second-rank spacelike tensors and the orthogonally projected component of 
vectors. For example, D^u^ = D((,u ) — (D c u c )h a b/3 and E/ a \ = h a h Eb - see Eq. Q. Also, overdots represent 
proper-time derivatives and primes conformal-time ones - see Eq. JHJ. 

3 The totally antisymmetric tensor e a b c , with e a b c u c = 0, is the Levi-Civita symbol of the 3-D space. 



2 



which are supplemented by the constraints 

D a E a = 11- 2u) a B a and B a B a = 2u a E a . (6) 

Here J a is the electric 3-current, \x is the electric charge and curLB a = e a bcD b B c by definition 
(with an analogous expression for curLE a ). The 3-current is related to the electric field via Ohm's 
law. For a single charged fluid the latter takes the form 

Ja = <>E a , (7) 

with q representing the electrical conductivity of the medium [7J. At the ideal magnetohydro dy- 
namic (MHD) limit the conductivity is very high, with <r — > oo. Then, the electric field vanishes 
and the currents keep the magnetic field frozen-in with the fluid. At the opposite end, where 
q — > 0, the currents are zero despite the presence of a finite electric field. Using the set ([Si)-© we 
can follow the evolution of electromagnetic fields in a variety of cosmological environments and 
study the implications of large-scale -B-fields, in particular, for the formation and the evolution 
of the structure that we see in our universe today. 

3 Primordial magnetic fields 

The detection of coherent magnetic fields in remote astrophysical systems increases the possibil- 
ity of significant .B-fields of cosmological origin. Early magnetogenesis, however, faces problems 
with both the strength and the size of the primordial field. Typically, magnetic fields gener- 
ated between inflation and recombination have too small coherence lengths and cannot seed the 
galactic dynamo. The reason is causality, which limits the size of the seed to the horizon scale 
at the time of magnetogenesis. This is typically much smaller than the 10 Kpc length required 
by the dynamo. A mechanism known as 'inverse cascade' can increase the correlation length by 
transferring magnetic energy to larger scales, but needs large-amounts of helicity [8jQ 

Inflation has long been seen as a solution to the scale problem, since it naturally creates super- 
horizon correlations. There is a serious strength problem however. Magnetic fields that survive 
an epoch of typical inflationary expansion are weaker than 10~ 50 G. This lies well below the 
dynamo requirements, which vary between 10 -12 and 10~ Gauss, in today's values, depending 
on the efficiency of the amplification and the cosmological model it operates in [21 [5]. This 
dramatic depletion is attributed to the 'adiabatic decay' of large-scale magnetic fields, which 
typically dilute as a~ 2 (a is the cosmological scale factor). To get an idea why, consider the 
expressions ([4]) and ([5]) of the previous section. These combine to give a wave equation for each 
component of the Maxwell field |9j. On a Friedmann- Robertson- Walker (FRW) background, the 
linear wave formula for the n-th magnetic mode reads as 

B'( n) + n 2 B {n) = -2KB {n) , (8) 

where Bi n ) = fl2 -^(n) is t ne rescaled magnetic vector, n is the eigenvalue of the mode, K = 
0, ±1 is the 3-curvature index of the FRW background and the primes indicate conformal-time 

4 See also Prof. Brandenburg's contribution to this issue. 
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derivativeso When K = the above reduces to the Minkowski-like expression 

B'( n) + n 2 B {n) =0, (9) 

which solves to give Br n \ oc a~ 2 at all times and on all scales. In practice, this result trans- 
lates into a current (comoving) strength below 1CP 50 G for typical 'inflationary' magnetic fields. 
Therefore, the adiabatic (B( n ) 00 a ~ 2 ) decay needs to slow down if we are to have inflation- 
generated -B-fields of astrophysical relevance. The effect is known as 'superadiabatic amplifica- 
tion' and it is usually achieved outside classical electrodynamics (e.g. see p~0j for a representative 
though incomplete list), unless FRW models with nonzero spatial geometry are employed 



4 Magnetic effects on structure formation 

The isotropy of the CMB strongly suggests that our universe was extremely smooth at recom- 
bination. Today, however, we see structure all around us and the obvious question is how this 
structure was formed. Gravitational instability seems to hold the answer, but we do not know 
the details yet. The ACDM model is the current 'concordance' scenario, which however excludes 
magnetic fields and has a rather large number of free parameters^ 

Studies of magnetised structure formation typically work within the ideal MHD approxima- 
tion and look at the effects of the magnetic Lorentz force on density inhomogeneities o These 
generally come in the form of scalar, vector and (trace-free) tensor distortions. The former are 
those commonly known as density perturbations and represent overdensities or underdensities in 
the matter distribution. Vector inhomogeneities are monitored by the curl of the density gradi- 
ent and describe rotational, vortex-like, density perturbations. Finally, tensor inhomogeneities 
correspond to shape distortions. Following [5], the scalar 

A = -D 2 p, (10) 

p 

describes linear density perturbations (D 2 = D a D a is the 3-D Laplacian operator). In a per- 
turbed, weakly magnetised, spatially flat FRW universe, the above evolves according to [5] 

A = 3wHA - (1 + w)Z + ^ c 2 (l + w)HB , (11) 

5 We use continuous harmonic eigenvalues, with n 2 > 0, in spatially flat and open FRW models and discrete 
ones, with n 2 > 3, when K = +1 g]. 

6 Friedmann models with non-Euclidean spatial geometry are only locally conformal to Minkowski space. On 
these backgrounds the adiabatic (B oc a~ 2 ) magnetic decay is guaranteed only on scales well inside the curvature 
radius. Near and beyond the curvature length, the magneto-geometrical term in the right-hand side of Eq. ((8} 
becomes important and can change the standard evolution of the B-field. 

7 The ACDM model assumes an accelerating universe dominated primarily by a cosmological constant (A) - 
or by dynamical dark energy acting as an effective cosmological constant - and secondarily by Cold Dark Matter 
(CDM). Baryons make only a small fraction (~ 5%) of the total matter. 

8 Beyond the linear regime, one should also account for microphysical processes and departures from the ideal 
MHD limit. As yet, however, studies of this nature are both sparse and patchy. 
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where Z = a 2 D 2 and B = (a 2 / B 2 )D 2 B 2 describe linear inhomogeneities in the expansion and 
the magnetic energy density respectively!^ To first order, these two variables propagate as 



11 c 2 1 

Z = -2HZ--pA + -c 2 (l + w)pB S —B 2 A - - c 2 B 2 B (12) 

2 F 4 aV 7P 1 + w 2 a v; 

and 

■ 4 ■ 4(c 2 - w)H . , , 

B = — r A + -H — A , (13) 

3(1 + w) 1 + w 

respectively. Also, w = p/p, H = a/ a is the background Hubble parameter, c 2 = p/p is the 
square of the adiabatic sound speed and c 2 = B 2 / p{l + w) is that of the Alfven speed. Finally, 
we note that B 2 <C p, given the relative weakness of the magnetic field. 

The system (llip -(fT3l) has analytical solutions in the radiation and the dust eras [5]. Before 
equipartition, when w = 1/3 = c 2 , H = l/2t, p = 3/ At 2 and c 2 = 3i? 2 /4p = constant, large-scale 
magnetised density perturbations obey the power-law solution 

A = C + CiH+f ^ + c 2 t~l c * + Crf 1 -! c * . (14) 

In the absence of the -B-field we recover the standard growing mode of A oc t, which means that 
the magnetic presence reduces the growth rate by 4c 2 /9. Well inside the horizon, the fc-mode of 
the density contrast oscillates like a magneto-sonic wave with 



A (fe ) oc sin 



l + 2 -^(^ J- 



(15) 



where = a/k is the perturbed scale and A# = 1/H the Hubble horizon [5]. Here, the magnetic 
pressure increases the effective sound speed and therefore the oscillation frequency. The former 
makes the Jeans length larger than in non-magnetised models. The latter brings the peaks 
of short-wavelength oscillations in the radiation density closer, leaving a potentially observable 
signature in the CMB spectrum |12j . 

When dust dominates, w = = c 2 , H = 2/3t, p = A/3t 2 and c 2 = B 2 /p oc t~ 2 / 3 . Then, on 
superhorizon scales, density perturbations evolve as [5] 

A = C + Ci^ 2 / 3 + C 2 t ai + C-zt a * , (16) 



with ai, 2 = - [1 ± 5^1 - (32/75)(c a X H /X k )l\/Q. In the absence of the 5-field we recover again 
the standard solution with ol\ = 2/3 and a 2 = — 1. As before, the magnetic presence slows 
down the growth rate of density perturbations. Also, the field's pressure leads to a magnetically 
induced Jeans length, below which density perturbations cannot grow. The magnetic Jeans 
scale, as a fraction of the Hubble radius, is [5] 

Aj ~ c a X H ■ (17) 



'Equation shows that B-fields are generic sources of linear density perturbations. Indeed, even when A 
and Z are zero initially, A will take nonzero values solely due to the magnetic presence. Note that only the 
pressure part of the Lorentz force contributes to the linear expressions (|lip and (|12[l . To account for the tension 
effects, one needs to allow for FRW backgrounds with non- Euclidean spatial geometry [5]. 
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Assuming a 5-field of approximately 10 -9 G, we find that Aj ~ 10 Kpc. Alternative, magnetic 
fields close to G, like those found in galaxies and galaxy clusters, give Aj ~ 1 Mpc. This 
lies intriguingly close to the size of a cluster of galaxies. Overall, the magnetic effect on density 
perturbations seems rather negative. Although l?-fields generate this type of distortions, they 
do not help them to grow. Instead, the magnetic presence either suppresses the growth rate 
of density perturbations or increases the effective Jeans length and therefore the domain where 
these inhomogeneities cannot grow. 

Magnetic fields also induce and affect rotational, vortex-like, density inhomogeneities. These 
are described by the vector W a = —(a 2 /2p)e a i, c D b D c p, which on an FRW background and after 
matter-radiation equality evolves as 

W a = -AHW a -\pWa + \ C a D 2 W a , (18) 

to linear order [5]. Defining A a = c a A# as the 'Alfven horizon', we may write the associated 
solution in the form [5] 

W {k) =Cit Ql +C 2 t a \ (19) 



with a± } 2 = — [5± J 1 — (48/9)(A a /Afc)g]/6. On scales far exceeding the Alfven horizon, A a <C A^ 

and the perturbed mode decays as W a t~ 2 / 3 . This rate is considerably slower than W oc t~^, 
the decay rate associated with magnetic- free dust cosmologies. Well inside A a , on the other 
hand, magnetised vortices oscillate like Alfven waves, with [5] 



W(jfc) oc t~ 5/6 cos 



2V3 /Aa 
9 \XkJo 



hit 



(20) 



Thus, the effect of the 5-field on a given vortex mode is to reduce its standard depletion rate. 
Analogous is the magnetic effect on u a , the vorticity proper. Overall, magnetised cosmologies 
rotate faster than their magnetic-free counterparts. In contrast to density perturbations, the 
field seems to favour the presence of vorticity. This qualitative difference should probably be 
attributed to the fact that the tension part of the Lorentz force also contributes to Eq. (|18p . 



5 Discussion 

In addition to scalar and vector perturbations, magnetic fields also generate and affect tensor- 
type inhomogeneities that describe shape-distortions in the density distribution [5j. An initially 
spherically symmetric inhomogeneity, for example, will change shape due to the magnetically 
induced anisotropy. All these are the effects of the Lorentz force. Even when the latter is removed 
from the system, however, the B-Held remains active. Due to its energy density and anisotropic 
nature, for example, magnetism affects both the local and the long-range gravitational field. 
The anisotropic magnetic pressure, in particular, leads to shear distortions and subsequently to 
gravitational-wave production. Overall, magnetic fields are a very versatile source. They are 
also rather unique in nature, since -B-fields are the only known vector source of energy. An 
additional unique magnetic feature, which remains relatively unexplored, is its tension. When 
we add to all these the widespread presence of magnetic fields, it makes sense to say that no 
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realistic structure formation scenario should a priori exclude them. It was probably thoughts 
like this that motivated some researchers to start adding magnetic fields into their numerical 
codes [13] , Hopefully their numbers will increase and we will soon have structure formation 
models with fewer free parameters and more physics. 
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